Lorentz invariant photon number density 
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A Lorentz invariant positive definite expression for photon number density is derived as the 
absolute square of the invariant scalar product of a polarization sensitive position eigenvector and 
the photon wave function. It is found that this scalar product is independent of the form chosen for 
the wave function and that the normalized positive frequency vector potential-electric field pair is 
a convenient choice of wave function in the presence of matter. The number amplitude describing 
a localized state is a 5-function at the instant at which localization and detection are seen as 
simultaneous. 



I. INTRODUCTION 

The concept of photon number density arises in the 
interpretation of experiments such as photon counting 
and the creation of correlated photon pairs in a nonlin- 
ear material. In spite of its relevance to experiment and 
the foundations of quantum mechanics, no relativistically 
invariant positive definite expression for photon number 
density exists at present. Mandel [l[ defined a coarse- 
grained number operator to count photons in a region 
large in comparison with the wavelength of any occupied 
mode of the field. However photons can be counted, and 
hence localized, in a photodetector much smaller than 
this wave length. Here I will show that the absolute 
square of the Lorentz invariant scalar product of a local- 
ized state and the photon wave function is a local positive 
definite number density. 

In nonrelativistic quantum mechanics the 1-particle 
number density is the absolute square of the real space 
wave function. The wave function can be obtained by 
projection of the state vector onto a basis of simultaneous 
eigenvectors of the position and spin operators. Although 
this was long thought to be impossible for the photon, 
we recently constructed a photon position operator with 
commuting components 0, H| and derived a photon wave 
function in this way HQ. The Landau Peierls (LP) wave 
function [|| and the positive frequency vector potential- 
electric field (AE) wave function pair were considered 
and a scalar product was defined to complete the Hilbcrt 
space. It was proved that the scalar product is invariant 
under the similarity transformation relating the AE and 
LP forms of its integrand, implying that they are equiv- 
alent when used for calculation of transition amplitudes 
and expectation values. 

Since position is an observable, the number ampli- 
tude to detect a particle at position r with spin cr 
equals the scalar product of the corresponding position- 
spin eigenvector with the particle's wave function. 
In real space this number amplitude, (r, (<)) = 
/ cPr'ipra- ( r ') ( r ';*)j i s equal to the wave function, 
^ct (r,t), only if the position eigenvectors are localized 
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states of the form -0 r . CT (r') = S 3 (r' — r) . While this is 
true for a nonrelativistic massive particle, it is not the 
case for the photon. A transverse vector is not of the 
simple (5-function form, and thus it is not the amplitude 
to detect a photon at r. In addition, photons are of- 
ten most conveniently described in terms of potentials or 
fields because of their simple relationship to the matter 
current density. The Fourier components of the vector 
potential and electric field describing a localized state go 
as u^ 1 ^ 2 where k is the wave vector and tOk is the angular 
frequency. They are not ^-functions even if their vector 
properties are ignored. 



Since the boost operator just generates a change of 
point of view that does not change the results of possible 
experiments 7], the scalar product that predicts these 
results should be a relativistic invariant. Newton and 
Wigner (NW) defined an invariant scalar product and a 
position eigenvector at the origin that is invariant un- 
der rotation and inversion Q. This and the displaced 
single particle states generated in k-space by the spatial 
translations exp (— ik • r) are orthonormal. Philips de- 
fined a Lorentz-invariant localized wave function but lost 
the orthonormality condition 9]. The NW wave func- 
tions themselves are not local in real space even in the 
simplest spin-zero case since they go as w^/ 2 - However, 
as argued in the paragraph above, the number amplitude 
is not in general equal to the real space wave function. 



In this paper, the relationship between the photon 
wave function and number density will be examined. In 
Section II, our previous work on the photon position op- 
erator and wave function will be summarized. In Section 
III a Lorentz invariant scalar product will be defined. 
The number amplitude for arbitrary polarization, equal 
to the scalar product of a position eigenvector and the 
wave function, will be found in Section IV. In Section V, 
interaction with matter will be considered, and photon 
number density will be discussed in relation to the recent 
and historical literature. 
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II. POSITION OPERATOR AND WAVE 
FUNCTION 

The photon position operator with commuting compo- 
nents and transverse eigenvectors can be written as [3| 



as 



?( Q >x) — 



(1) 



where V is the k-space gradient and D = 
cxp (— iSkX) ex P {—iSz4>) ex P is the rotation 

matrix with Euler angles </>, 9, x- The operator D rotates 
the unit vectors x, y and z into the spherical polar 
unit vectors <f> and k and then rotates 9 and <p about 
k by x(k). Thus r^ a,x ^ rotates the transverse and 
longitudinal unit vectors to the fixed directions x, y 
and z, eliminates the factor (cj^) , operates with iV to 
extract the position information from the phase of the 
wave function, and then reinserts the factor (cok) and 
returns the unit vectors to their original orientations. 
The eigenvectors of (HJ form a basis for the LP wave 
function if a = 0, and for a wave function proportional 
to the vector potential if a = —1/2 and its conjugate 
momentum if a = 1/2. 

In the Schrodinger picture (SP) the simultaneous 
eigenvectors of the position and helicity operators, 



Vfe x) (k) = (a,*) 



a exp (— zk • r) 



(27T) 



.3/2 



(2) 



form a basis for the Hilbert space, TL. The transverse 
unit vectors are 

4*> (k) = i= (0 + iat) exp (-ia X ) (3) 

for circular polarizations a = ±1 and the longitudinal 
unit vector is eo = k. The real linear polarization unit 
vectors 



e W - _ U ix) + e ix) \ 
1 ~ v/2 V ' 



(4) 



_L( e (x)_ e (x)> 



,{x) 



give eigenvectors of r'"^) no t the helicity operator. 
While all x arc needed to interpret an experiment that 
measures polarization, only one basis is needed to de- 
scribe the photon state, and the x = definite helicity 
basis will be used here for simplicity. In this basis, the 
probability amplitude for the state with polarization 
incorporates the phase factor exp (— icrx)- 

The wave function was obtained in [4| as the projection 
of the quantum electrodynamic (QED) state vector onto 
a basis of position- helicity eigenvectors. If the operator 
aj. (k) creates a photon with wave vector k and circular 

polarization ei , the 1-photon k-space basis states are 
|k, a) — a\ (k) |0) where |0) is the vacuum state. Cre- 
ation operators for a photon at position r can be defined 



(«)t 



i , a exp(-zk-r) 

(wfc) 575 — e. 

(2tt) 3/2 



(o)* 



(k)4(k). (5) 



Any state vector can be expanded in Fock space as 
I^W) — S^Lo Cn M) wne re the 1-photon term 
l*i(*)> = E«r/ d 3 fc (k,cr|* (<)) \k,a) is completely de- 
scribed by the probability amplitude 



(k, cr|* (t)) = C(j (k) exp (-iujkt) . 



(0) 



The projection onto the definite helicity position eigen- 



$2 



* (t) ) , is a six compo- 



vectors, (r, t) = ^0 

nent wave function whose dynamics is described by a di- 
agonal Hamiltonian 0, @j. The expectation value can 
be evaluated using [a a (k) , tv (k')] = <$ CTi(T /<5 3 (k — k'). 
Here the focus is on the scalar product and polarization 
will be summed over to give the 3-vector wave function 
* (q) (r, t) where 



*W (r,t)= /" d 3 fcc CT (k) 



(7) 



xe^(k)(w fc ) — ^ , 



(2*)° 



*W(r,t)= J] *£ a) (',*)• 



The k-space wave function will be defined as 

*W(k) =(V (k)e< )(k)Kr, (8) 
* (q) (k) = ^ *i Q) (k) • 



er=±l 



III. INVARIANT SCALAR PRODUCT 

To prove the invariance of the scalar product that com- 
pletes Tt, 4- vector notation and the Lorentz gauge will 
be used. The contravariant energy-momentum 4-vector 
is hk where k = k v = (ujk / c, k) and the covariant 4- 
vector is k^ = g^ v k v where g is a diagonal tensor with 
goo = — 1 and gu = 1 for i = 1 to 3. In addition to the 
transverse and longitudinal eigenvectors, scalar position 
eigenvectors can be define d I l0l| and the 4-potential can 
be used as wave function [ll| . In the Lorentz gauge the 
momentum conjugate to A(.+)f is FEW" = eodA<-+^/dt 
where the superscript (+) denotes the positive fre- 
quency part. This is equivalent to ip^ 1 / 2 ^ — idip^ 1 ^ 2 " 1 jdt 
that follows from ([7]). Since only positive frequencies 
arise in the wave function, the invariant volume integral 
can be written as 



d 4 kd (Jfe"jfe„ + k 2 ) 6 (oj k ) = c 



d 3 k 
2^ 



(9) 



where the Heaviside step function (uJk) is invariant un- 
der proper Lorentz transformations and k = mc/h = 
for a free photon. 
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The real space wave function satisfying (x) — 

V / W^A(+) (x) [| is a 4-vector. Use of © in Q for 
a = —1/2 gives 

(x) = Aj!fcexp(to)^ (1/2) 
./ (2tt) 3/2 

for physical states with exp (— iwfct) time dependence. 

This implies that 'J^ 2 '' (fc) must be a 4-vector. Thus 
the scalar product of |$i) and , 

<*l|*l> = / — $(V2)M* (fc) ^(1/2) (fc) ? (n) 

is an invariant. This is the form normally used in 
field theory, incorporating the metric w^T 1 . Alter- 
natively the integrand $( 1 / 2 )' 1 * (fc) *p 1/2) (fc) /w fe can 
be written as the product of LP wave functions, 
or as the potential-conjugate momen- 
tum product, $(-!/2)m* (fc) vl/^ 1 / 2 ' to give the invari- 
ant scalar product (jTTJ) as 



(12) 



In the Lorentz gauge there are longitudinal and scalar 
photons in addition to the observable transverse photons. 
Scalar and longitudinal photons can be dealt with by in- 
troducing an indefinite metric fl2| . However the Lorentz 
gauge was only needed to prove invariance and its further 
use here is not required. 

In any specific reference frame it is possible to make a 
gauge transformation to the transverse (Coulomb) gauge. 
Eq. (fl"2"|) then reduces to 



= J d 3 k&- a > (k) • *W 



(k). (13) 



The conjugate momentum is n' + )= — eoEj_ for the 
minimal coupling Hamiltonian (_L denotes the transverse 
part), so the a — ±1/2 conjugate pair are the vector 
potential and the electric field (AE). Since the {u>k) ±a 
factors cancel and the dot product of like unit vectors is 
unity Eq. (|13p can also be written as 



f d 3 kd:(k)c„(k) 

Li J 



(14) 



<T=±1 



where the d's are the coefficients in the expansion of |$i) . 
This is just the QED scalar product. 

For any operator, O, describing an observable, 
O^a (k) must be in 7i and the eigenvalues of O must be 
real. This is obviously the case for the momentum and 
energy operators, ?ik and hw^. For the position operator, 
using ©, 

yKx)^*) (k ) = (wfc )« e (x) ( k ) iV fa (k, t) exp (ia X )} 

(15) 



since {wk) a (k) is an eigenvector of r( Q,x ) with eigen- 
value zero. Thus (fl~5)) is in Ti. Using the scalar product 
([TT]) the wave function is just vl^ 1 / 2 )^ (fc) and it can be 
proved as in [l3| that rU/ 2 '*) is Hermitian. If (fT3"|) is 
used, it can be proved using CE3 that (r(- a ^)$i|*!) = 
($ 1 |r(«.x)^ 1 ^* which impl ies that r'°' x ' is Hermitian and 
must have real eigenvalues and orthonormal eigenvectors. 
This follows the usual rules of operator algebra since D 
is unitary and iV is Hermitian. The a = ±1/2 basis is 
related to the a = basis by r( Q,x ' = (usk) a r^ '*) (u^.) - ". 
This similarity transformation preserves the scalar prod- 
uct and hence expectation values and the reality of the 
position eigenvectors [l,[l4j]. The wave functions 1 
and >I'ct 1 ' /2 ' ) form a biorthonormal pair and the position 
operator satisfies r( -1 / 2 ' x )t = fw 2 '*), so it can be called 
pseudo-Hermitian [l4l |. 

In real space, the most convenient form for the scalar 
product is not obvious. If dTTJ) is transformed directly, 
the metric factor uj^ 1 must be replaced by the inverse of 
the Hamiltonian operator, implying a nonlocal integrand 
6] . Eq. (fT3")) transformed to real space is 

($1!*!) = J Sr&- a > (r, t) • #W (r, t) , (16) 

as can be verified by substitution of {7]) and integration 
over d 3 r to give <5 3 (k — k') and then (JT3J) - The integrand 
of (|16[) is local, making it a useful form of the scalar 
product. 

It can be seen by inspection of Eqs. (fTTj) to (fT6|) that 
the scalar product is unaffected by the change of metric 
from (fTTj) to (|12[) . and is invariant under the similarity 
transformations between the AE and LP forms of the 
wave function and under the unitary transformation be- 
tween r-space and k-space. It is also invariant under the 
unitary transformations exp (— iSkx) to x ¥^ bases. The 
AE form of the wave function is preferable in most ap- 
plications, since the relationship of the LP wave function 
to matter source terms is nonlocal [l5j], but the choice is 
a matter of convenience. 



IV. PHOTON NUMBER AMPLITUDE 

The probability amplitude for a measured result is 
the amplitude for a transition from the photon state 
described by c a (k) to a final state that is an eigen- 
vector of the operators representing the experiment. 
For a measurement of momentum fik and polariza- 
tion this final state is described by d a > (k') = 
b a ' ,o exp (—iax) S 3 (k' — k). Substitution in (fl"4|) gives 
the probability amplitude c a (k) exp {iax) where the fac- 
tor exp {iax) rotates the polarization about k by — X- 
The helicity is an invariant and a Lorentz transfor- 
mation just changes ^ @, so the probability to de- 
tect a photon with definite helicity and wave vector 
is invariant. For a measurement of the linear polar- 
ization e^, given by ((4]), the momentum eigenvectors 



4 



are d a (k') = exp (-icrx) S 3 (k' - k) /y/2 and d a (k') = 
i(-l) cr exp(-icrx)<5 3 (k'-k)/V2 for j = 1 and 2 re- 
spectively, so (fT4|) gives the probability amplitude to de- 
tect a photon with wave vector k and this polarization 
direction as 

cri (k, x) = [ci (k) exp (ix) + c_i (k) exp (-«%)] , 

(17) 

cm (k, x) = -j= [ci (k) exp (ix) - c_i (k) exp (-ix)] • 

If x = the measured polarization directions are 9 and 
0, and cri and cr2 are the amplitudes for transverse 
magnetic and transverse electric modes respectively fl6| . 
Probability density is the absolute valued squared, so the 
probability to detect a photon with definite helicity is \~ 
independent, but a linear polarization measurement is 
sensitive to phase. For c a oc exp (— icrx'), cri oc cos Ax 
and CL2 oc sin Ax where Ax = x' ~ X 1S the polarization 
angle of the photon relative to the polarization measured 
by the apparatus. 

The physical states with circular polarization lo- 
calized at position r at fixed time t have the k-space 
amplitudes 

, , \ c exp (-icrx - «k ■ r + iu k t) 
d <T >(k,T,t,x)=5 <T >, <7 /2 . (18) 



(2tt 



This is the probability amplitude for a position eigenvec- 
tor if the phase x — > X — crcjfct in |T]). The amplitude for 
a photon in state c a (k) to make a transition to this state 
is given by substitution of (fT81 ) into (fT4")l as 



Cff(M) = /d » fc 5E(iox + &-r-fc*i) 



(27T) 



3/2 



Cr(k). (19) 



It equals the inverse Fourier transform of 
c CT (k) exp (icrx - K^i). 

As a consistency check, consider a photon with po- 
larization localized at r' at time t' . The probabil- 
ity amplitude for this state is c a (k) = d a (k, r',t',x) 
and (|19p can be integrated to give the number amplitude 
c CT (r, t) — <5 CTi(T '(5 3 (r — r') if t — t', that is if localization 
and measurement are seen as simultaneous. All k's are 
included with equal weight in (fT5)) and exact localiza- 
tion, which is limited by the Hegerfeldt theorem [17|, is 
only possible because of interference effects between the 
converging and the diverging wave. Eq. (|18p describes 
a situation where there is no physical detector at r, so 
absorption is followed immediately by re-emission. 



V. CONCLUSION 

The most widely accepted photon wave function in the 
current literature is the Glauber photodetection ampli- 
tude [3, e|l +) (r,t) = (o E(r,t) where E(r,i) is 



the Heisenberg picture (HP) electric field operator and 
Ej^' is proportional to SI^ 1 / 2 ) given by (7}. This can 
be combined with <I>( -1 / 2 ) to give a scalar product with 
a local integrand and a complete description of photon 
quantum mechanics. Since the scalar product is invari- 
ant under unitary and similarity transformations and 
the prediction of the experimental results requires only 
the scalar product, the choice of wave function can be 
based on convenience. If desired, the remaining fields 
DW = e E(+) +P(+), H(+) and B« = Mo H(+) +M« 
can defined similarly 0, E^] where |0) is the matter- 
field ground state and P and M create and destroy 
matter excitations. Since the HP field operators satisfy 
classical dynamical equations, the equations satisfied by 
these positive frequency fields are identical in form to 
Maxwell's equations and describe the dynamics of the 
l-polariton state. 

Evaluation of the photon number amplitude (|19[) is 
straightforward since it just requires integration of the 
scalar QED probability amplitude. Polarization unit vec- 
tors as a function of k can be selected for convenience, 
for example, x = — 4> gives e^~^ (fez) = (x + icry) /y/2 
for a paraxial beam propagating parallel to z. Number 
amplitude itself has intrinsic physical significance. It is 
number amplitude that, after Schmidt diagonalization, 
was used as a photon wave function by Chan, Law and 
Eberly [2l[. Also, number amplitude is proportional to 

Ej^ to a good approximation. The positive frequency 
electric field of a localized photon state drops off as r~' 7 1 2 
where r is the distance from the point of photon local- 
ization (20| . Thus photon number density, equal to the 
absolute square of (|19p , is an excellent approximation to 
the Glauber photodetection probability, proportional to 

Number density may prove to be of fundamental im- 
portance. Cook pointed out that photon current density 
cannot be made precise within the Glauber theory, so 
that it fails to provide a complete description of pho- 
ton transport [15| . Calculation of the probability density 
to detect a photon as the absolute square of the scalar 
product of the 1-photon wave function and a position 
eigenvector using (|19p follows the usual rules of quan- 
tum mechanics. There is no such connection between 
Glauber photodetection theory and the basic rules of 
quantum mechanics. It is an ongoing goal of this au- 
thor to thoroughly understand the relationship of pho- 
ton number density to photon fields, energy density, the 
interpretation of photon counting experiments, and fun- 
damental issues relating to photon localizability. 

In summary, photon number amplitude in real space 
was calculated as the invariant scalar product of a lo- 
calized state with definite polarization and the 1-photon 
wave function. It gives an excellent approximation to the 
Glauber photodetection probability, with the advantage 
that the orthonormal position eigenvectors lead to mu- 
tually exclusive probability densities and an integrated 
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probability of unity. When described in terms of this 
number amplitude, the localized basis states combine 
the orthonormality of the NW states with the Lorentz 
invariance of Philips' localized states [j| Q. The posi- 
tive frequency vector potential and electric field provide 
a complete description of photon quantum mechanics, 



including interaction with matter. This should make it 
possible to test the relevance of the photon number con- 
cept to experiment. 
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